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ABSTRACT 

The Delaware Method of rating baffled shell and 
tube heat exchangers is used as the basis for a simplified 
mathematical model for predicting shell side heat transfer 
coefficients and pressure drops. Multiple non-linear 
regression analysis is used to express the heat transfer 
coefficient and pressure drop as a product of terms 
involving the five basic geometric parameters of a shell- 
and-tube heat exchanger, over a suitable range of these 
variables . 

A brief sensitivity analysis of heat transfer per 
unit pressure drop is presented, along with a proposed 
routine for determining exchanger size when the shell- 
side coefficient is the controlling factor. A comparison 
is made between this method and the Delaware Method. 
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NOMENCLATURE 

? 

heat transfer area for tube side, (ft ) 

2 

A heat transfer area for shell side, (ft ) 

0 

c s specific heat of shell side fluid, (Btu/lbm - °F) 

shell inside diameter, (in) 

D _ . diameter of outer limit of tube bundle, (in) 
ot£ ’ J 

d Q tube outside diameter, (in) 

d^ tube inside diameter, (in) 

F fraction of tubes that are in crossflow 

c 

F , fraction of total crossflow area, available for 

s P bypass flow around tube bundle 

f. gravitational friction factor for flow across 

1 an ideal tube bank 

O 

g gravitational conversion constant, 4.17 x 10 

c Ibm-ft/lbf-hr* 

h^ tube side heat transfer coefficient (Btu/hr- ft^ - °F) 

ideal shell side heat transfer coefficient, 
ldeal (Btu/hr-ft z -°F) 

h exchanger shell side heat transfer coefficient, 

0 (Btu/hr-ft 2 -°F) 
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J, shell side heat transfer coefficient correction 

factor for bundle bypass effects 

J shell side heat transfer coefficient correction 

factor for baffle configuration effects 

J. shell side heat transfer coefficient correction 

^ factor for baffle leakage effects 

J shell side heat transfer coefficient correction 

factor for adverse temperature gradients 

J shell side heat transfer coefficient correction 

factor for the effect of unequal baffle spacing 

Colburn j -factor for ideal tube bank 

K, thermal conductivity of tube side fluid, 

D (Btu/hr-ft-°F) 

K thermal conductivity of shell side fluid, ' 

s (Btu/hr- ft-°F) 

K w thermal conductivity of tube wall, (Btu/hr- ft -°F) 

•^£m logarithmic mean temperature difference 

L effective tube length between tube sheets, (ft) 

l baffle cut distance from baffle to shell inside 

diameter, (in) 

£ s baffle spacing, (in) 



baffle spacing at inlet- and outlet of exchanger 
respectively 

number of baffles in exchanger 

number of tube rows crossed during flow through 
one crossflow section 

number of effective crossflow rows in each 
window section 

total number of tubes in exchanger 

pressure drop for flow across one ideal crossfl 
section, (lbf/in 2 ) 

pressure drop for flow through one ideal window 
section, (lbf/in 2 ) 

tube pitch, (in) 

tube pitch normal to flow, (in) 

tube pitch parallel to flow, (in) 

total heat transferred in heat exchanger, 



(Btu/hr) 



pressure drop correction factor for effect of 
bundle bypass 



pressure drop correction factor for effect of 
baffle leakage 



pressure drop correction factor for effect of 
unequal inlet and outlet baffle spacing 
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m 



sb 



tb 



w 



wg 



wt 



U 



W 



w. 



sb 



tb 



st£ 



shell side Reynolds number 

crossflow area at centerline for one crossflow 
section, (in 2 ) 

2 

shell -to-baffle leakage area for one baffle, (in ) 

2 

tube- to-baffle leakage area for one baffle, (in ) 

2 

area for flow through window, (in ) 

2 

total baffle window area, (in ) 

2 

window area occupied by tubes, (in ) 

overall heat transfer coefficient base on shell 
side heat transfer area, (Btu/hr-ft 2 -°F) 

mass flow rate of shell side fluid, (lbm/hr) 

mass flow rate of tube side fluid, (lbm/hr) 

diametral clearance between shell and baffle, 

(in) 

diametral clearance between tube and baffle, 

(in) 

approximate diametral clearance between shell 
and outer tube limit, (in) 

shell side fluid viscosity at bulk temperature, 
(lbm/hr-ft) 



shell side fluid viscosity at surface temperature, 
(lbm/hr-ft ) 

shell side fluid density, (lbm/ft^) 
baffle cut angle, (rad) 
shell side Prandtl number 

tube side fluid viscosity at bulk temperature, 
(lbm/hr-ft) 



1 

2 

3 
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I. INTRODUCTION 

» 

A. BACKGROUND 

The aim 'of any heat exchanger rating method is to 
predict, as accurately as possible, the performance of a 
heat exchanger with respect to total heat transferred and 
pressure drops encountered on both the shell and tube sides. 

The Delaware Method for Rating Baffled Shell and 
Tube Heat Exchangers (1) was developed from its inception 
in the late 1940’s through the early 1960’s as a result of 
numerous experiments conducted at the University of 
Delaware. This method is applicable to many variations 
on the baffled exchanger: multiple passes, finned tubes, 

various tube orientations, a wide range of Reynolds 
numbers, etc. As Bell has stated (1), the method does 
have its uncertainties, as all methods do. However, this 
method is widely used and easy to understand. 

Tinker (2) and Devore (3) propose similar methods 
for rating baffled exchangers, but, in the opinion of 
this author, are not as straightforward as the Delaware 
Method. Bell himself has proposed routines for sizing 
baffled exchangers based on the Delaware Method 
(references 1 and 11) , and this method has been used 
extensively in computer programs for "optimal” exchanger 
determination. References (9) and (10) are two such 
examples . 
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The basic means o£ analysis in the Delaware Method 
takes into consideration five fluid streams and three 
heat transfer zones in the baffled exchanger, as illustrated 
in Figures 1 and 2. As can be seen by these simplified 
illustrations, the flow pattern for shell side fluid is 
complex at best. Reference (1) contains the complete set 
of equations for relating the heat transfer coefficient 
and pressure drops to the geometry of the exchanger. This 
reference also utilizes extensively graphs and tables 
for determination of some parameters. Mueller (7) has 
developed equations for many of these graphical relations. 
Appendix A is a compilation of the basic . equations of the 
Delaware Method and those derived by Mueller. This 
compilation does have its limitations as noted in Section I.C. 




FIGURE 1 Schematic Diagram of Streams on 
Shell Side of Baffled Heat 
Exchangers . 
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FIGURE 2(a) Region o£ Crossflow Between Baffle 
Tips in. the Central Baffle Spacing 




FIGURE 2(b) Flow Region Considered for Window Flow 




FIGURE 2(c) Flow Region for the End Baffle Spacings 
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B. PURPOSE AND METHODOLOGY 

The purpose of this study is to provide an insight 

into the relations between shell-side heat transfer 

coefficient and pressure drop and the five basic 

exchanger parameters: Tube diameter -- d Q ; pitch-to- 

diameter ratio -- P/d : baffle cut ratio -- £ /D. ; baffle 

o C l 

spacing ratio -- £ s /D^; and shell diameter -- D^. 

It is anticipated that a simplified model for 
determining the heat transfer coefficient will be of 
considerable use to the designer as a "first cut" in 
sizing an exchanger for a given job. Admittedly, this 
cannot replace the sophisticated computer programs 
available for the determination of "optimal" exchanger 
configuration for a given application. However, the 
method developed here will be a step in the right 
direction for initial selection of parameters and will 
provide some "feel” for the effect of geometry variations 
on overall heat exchanger performance. 

This analysis was conducted utilizing the basic 
equations of the Delaware Method as modified by the 
assumptions and approximations enumerated in Section I.C. 
Utilizing these equations, those factors affecting heat 
transfer and pressure drop which are a function of 
exchanger geometry were calculated over a suitable range 
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of variables. The resulting values were correlated 
utilizing four FORTRAN subroutines (b) : CORRE, ORDER, 

MINV and MULTR. The resulting non-linear correlation 
yields these factors expressed as a product of the 
five basic exchanger parameters, each with its own 
regression coefficient as an exponent, and a single 
correlation coefficient. These regressions provided 
acceptable correlation results, the maximum overall 
error being on the order of five percent. 

C. LIMITATIONS/ASSUMPTIONS/APPROXIMATIONS 

1. Limitations 

In this analysis, no attempt has been made to 
correlate data for the shell side heat transfer co- 
efficient where sealing strips are used to limit the 
bundle bypass flow. Nor has an attempt been made with 
regard to enhanced surfaces for heat transfer, i.e. 
finned tubes. Finally, only the single pass condition 
for the tube side fluid is taken into consideration. 
Admittedly, these restrictions limit the scope of 
this study with regard to overall applicability to 
heat exchanger design. However, it is hoped that this 
analysis will provide insight into the basic relations 
between shell side heat transfer coefficient and the 
fundamental exchanger goemetric parameters. 



-17- 



2. Assumptions 

It is realized that manufacturing considerations 
effect significantly the design, construction and final 
form of any heat exchanger. Reference (5) sets standards 
for the tube- to-baffle and baffle- to-shell diametral 
clearances for various applications. These values are 
predominantly of the "step function" type, in that they 
can not be expressed as a continuous dependent variable 
based on any of the independent exchanger parameters. 

For the purposes of this study, both clearances were 
assumed constant throughout the range of parameters 
considered. The values used for these clearances are: 



Shell- to-baffle clearance: 6 ^ 

Tube- to-baffle clearance: 5^.^ 



0.175 inches 
0.03125 inches 



These assumptions again allow concentration on the effect 
of major exchanger parameters. 

The second major assumption is the neglect of 
unequal baffle spacing for inlet and outlet passages. 

This assumption is made in both heat transfer coefficient 
and pressure drop calculations. Unequal baffle spacing, 
in an exchanger where the ratios 



l 



s 






l 



s 



,o 



l 



s 



£ 



s 



1.5, 
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contributes to an approximate 2-3% reduction in the over- 
all heat transfer coefficient for a sufficiently large 
number of baffles. Where pressure drop is concerned, 
the assumption must be made with the realization that 
the end zones be considered as crossflow areas identical 
to the zones between the baffles. Equations (19) and 
(25) of Appendix A reflect these changes. 

3. Approximations 

In order to implement any numerical analysis of 
the type used in this study, it is necessary to derive 
relations from data provided by the manufacturing 
community and from curves of friction factor and 
Colburn j -factor. 

For the former, references (1), (5) and (8) were 
used for two such derivations. The number of tubes in a 
heat exchanger of particular shell diameter was corre- 
lated as a function of shell inside diameter, pitch- 
to-tube diameter, and tube diameter, for the single- 
pass condition. These results are: 

For square/rotated square tube layouts: 



N 



t 



0 17^fi 





2 



0.5559 



( 1 ) 
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For triangular tube layouts: 



N 



t 



0.40759 




2.1797 





( 2 ) 



These approximations give results within 51 over the range 
of shell diameters from 10 to 60 inches. The clearance 
for bundle bypass flow (D^ - D Qt ^) was correlated as a 
function of shell diameter also, resulting in the 
approximation : 



! . = (D. - D .„) = 0.5817 D. 

stJt i otl J 1 



0.5066 



(3) 



Figures 3, 4, and 5 are graphs of friction factor 
and Colburn j -factor versus shell-side Reynolds number 
for the various tube orientations at varying pitch-to- 
tube diameter ratios. Correlation results are presented 
in Tables 1 and 2 for friction factor and Colburn j -factor, 
respectively, for the separate ranges of Reynolds numbers 
shown. 

Table 3 presents relations between tube pitch, 
tube pitch normal to fluid flow and tube pitch parallel 
to fluid flow for the various orientations. These 
relations are used in equations (1) , (3) , and (4) of 
Appendix A in calculating the various factors associated 
with heat transfer and pressure drop calculations. 



- 20 - 




FIGURE 3 Ideal Tube^ Bank * j- 7 - af^ fr* Factors for 90° In-Line Layout 
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FIGURE 4 Ideal Tube Bank j. and f. Factors for 45° Staggered Layout 
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P = tube pitch 

P = tube pitch normal to flow 

Pp = tube pitch parallel to flow 



Tube 

Orientation 


P 

n 


P 

P 




square 


1.0 P 


1.0 


P 


rotated 

square 


/2/2 P 


/ 2/2 


P 


triangular 


P/2 


/3/2 


P 



TABLE 3 Relations Between Tube Pitch 
P, and Pjj, Pp for Various Tube 
Orientations 
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D. NEGLECTED FACTORS J r and J s 

For preliminary estimation and sizing as proposed 
by this simplification, the factors J r (correction due to 
adverse temperature gradients) and J s (correction due to 
unequal inlet and outlet baffle spacing) have been 
neglected. 

To account for the effect of J s , the final heat 
exchanger length must be known, which allows determination 
of the end baffle spacing. J s may then be calculated 
from: 



J s ■ 






(. 4 ) 



(• 4 ) 



S y 0 



(Nw_ 



b-1) 



+ £ 



s ,1 



+ £. 



s ,0 



where 



* 







/£. 



s , i' "s’ 



* 



£ 



s , o 




( 4 ) 



The initial value of h Q must then be multiplied by J g to 

give a more accurate estimate. 

The factor J is calculated only in laminar flow 
r 

(Re s < 100) from: 

1.51 

J r * 

(N c ) 



0.18 



( 5 ) 
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Should this condition (Re s < 100) exist, again adjust h Q 

by multiplying by J . Figures 6 and 7 show graphically 

* 

the relations for J vs l • and J vs N . These graphs 

^ ^ j JL X w 

may be used in lieu of actual calculations. 



- 28 - 




FIGURE 6 



Correction Factor for Unequal Baffle Spacing, 

J vs £* . ( = £* 

s s , 1 s , o' 
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FIGURE 7 Correction Factor for Adverse Temperature 

Gradients, J vs N 
’ r c 
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II. CORRECTION FACTORS INFLUENCING THE 
HEAT TRANSFER COEFFICIENT 



A. BUNDLE BYPASS FACTOR -- J b 

The bundle bypass factor, , accounts for the 
effect of fluid flow between the shell and outer tubes 
of the tube bundle. This portion of the flow bypasses 
the main crossflow section of the tube bank and has a 
significant effect on heat transfer. Equations A-4, 

A-5, and A-18 are used to calculate J^, where the shell- 
to- tube-bundle clearance (D^ - D has been approxi- 

mated by equation (3). Equation A-18 was developed by 
Mueller (7) where for laminar/transition flow C^ = 1*35 
and for turbulent flow C bh ~ 1 ' 25, For purposes of this 
analysis only two regimes of Reynolds numbers are 

considered: 1 < R P < 10* 5 for laminar/transition flow 

e s 

and 10^ < R e < 10$ for turbulent flow. 

The non-linear regression results for are 

presented in Table 4 and a representative curve is shown 

in Figure 8. As can be seen, can be represented as 

a function of shell diameter D- , tube diameter d , and 

1 o 

pitch- to-diameter ratio P/d . 

From Table 4, it can be readily seen that tube 
diameter has little effect on J^. Over the range of 
d Q considered (0 . 5<d Q <l . 5) , the change in due to 
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increasing (or decreasing) d Q from 1.0 to 1.5 (or to 0.5) 
is less than 1%. 

Figure 8 is applicable only for turbulent flow 
over square tube orientations. However, we may estimate 
for other flow/orientation conditions easily from 
Table 4 and can estimate changes in in the following 
fashion. If is desired for rotated- square layouts 
in turbulent flow, we may form a ratio: 



3. 3. 

\ _ CC 0 (P/d 0 ) 2 (Dp 3 ]r 

cl 3. 

J b s CC 0 (P/d 0 ) 2 (Dp Ps 




(P/d 0 ) 




( 6 ) 



s 



This simple relation yields the result that rotated- square 
layouts, for P/d Q = 1.5, = 30, result in a bundle 

bypass factor approximately 6% higher than square tube 
orientations . 

Such insight into the intricacies of the Delaware 
Method can allow the designer to make reasonable "fine- 
tuning" corrections to a preliminary design without the 
need for costly (and sometimes cumbersome) computer models. 
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J, = c 
b o 


a i 


(P/d Q ) a2 CD i ) 


a 3 




R 

e s 


Layout 


C 

0 


a i 


a 2 


3 -, 

O 




Square 


. 335 


- .0137 


. 5249 


.1806 


1 - 10 3 


Rotated- 

Square 


.416 


-.0119 


.4358 


.1493 




Triang . 


. 496 


- .0099 


. 3503 


. 1194 




Square 


. 369 


- .0125 


.4789 


. 1648 


10 3 - 10 5 


Rotated- 

Square 


. 445 


- . 0108 


. 3976 


. 1362 




Triang . 


. 527 


-.009 


. 3197 


.1089 



TABLE 4 Correlation Results for Baffle 
Configuration Factor, 







FIGURE 8 Bundle Bypass Factor J b vs P/d Q 
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B. BAFFLE CONFIGURATION FACTOR -- J c 

Equations A-2 and A-16 represent the baffle con- 
figuration factor, J c , which amounts for changes in heat 
transfer due to those tubes physically located in the 
baffle window and not subjected to full corssflow from 
the main stream of fluid. Equation A-16 is a linear 
approximation to the graph of J c developed by Bell. Using 
equation (3) for D Qt £> J c can be reasonably approximated 
as a product of shell diameter and baffle cut (or baffle 
cut ratio) . This result is independent of flow condition 
and tube orientation. Regression analysis results in 
the following equation for J c : 

-0.5467 -0.0184 

J c = 0.492 (Jl c /D i ) (D i ) (7) 

Figure 9 is a family of curves plotting i /D. vs J at 

C 1 c 

various values of D. . 

i 

As is easily seen from this figure, the dependence 

of J c on shell diameter is minimal when compared with 

that of the baffle cut ratio. For the range of baffle 

cut ratios considered in this analysis (0.2 < £ c /D^ < 0.4), 

the heat transfer coefficient may vary from a 15% increase 

to a 25% decrease. Bell (1) states that for "well-designed" 

exchangers, J c should be about 1.0, which corresponds to 

a value of £ /D - of about 0.25. 
c' 1 
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FIGURE 9 Baffle Configuration Factor J c vs Baffle 
Cut Ratio £ /D^ 
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It is interesting to note that in the optimization 
routines of references (9) and (10), the baffle cut ratios 
for ’’optimal" design were 0.25 and 0.2, respectively. 
Admittedly, the objective functions for these optimizations 
were slightly different (heat transfer area and total 
tube length, resp.), but the similarity of results is 
significant . 

C. BAFFLE LEAKAGE FACTOR -- J£ 

The baffle leakage factor, J^, is by far the most 
complex of all the factors affecting the overall heat 
transfer coefficient. It accounts for the effect of 
leakage through the baffle due to clearances between 
baffle and shell and between tube and baffle. In this 
regard, it is a function of all five basic exchanger 
parameters. Equations A-6, 7, 8, 15, 16, 17 must all 
be used in the calculation of J^, where some additional 
dependence is noted as a function of tube orientation. 

The non-linear regression results for are 
presented in Table 5. Numerous graphs and families of 
curves can be generated using these results. However, 
for the sake of simplicity, Figures 10a, b, and c are 
represenative of the effect on of several parameter 
variations. Bell states that the value of should 
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3. 3. 

J i ■ c 0 C<V 1 ( p/d c’ 2 


c/°/ 3 


3 A 

/d ± d 4 


^ c 

CD ± ) 5 




Layout 


C o a l 


a 2 


a 3 


a 4 


a 5 


Square 


0.572 0.1296 


0.7745 


0.0975 


0.2708 


0.0718 


Rotated- 

Square 


0.613 0.1177 


0.7136 


0.0856 


0.2339 


0.0643 


Triangular 


0.656 0.1127 


0.6679 


0.0767 


0.2051 


0.0516 


TABLE 5 


Correlation 

Factor, 


Results for 


Baffle 


Leakage 
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Square Tube Orientation 

d = 0.5, D . = 30 ; 
o ’ 1 




FIGURE 10a Baffle Leakage Factor J z vs P/d Q 
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FIGURE 10b Baffle Leakage Factor vs P/d Q 
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FIGURE 10c 



Baffle Leakage Factor J £ vs P/d Q 
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be 0.7 -- 0.8 for well-designed exchangers. From 
Figure 10, one can easily read an appropriate value 
of the necessary parameter that will be a reasonable 
starting point for heat exchanger design. For other 
shell diameters and tube diameters, the relations of 
Table 5 may be suitably manipulated to arrive at 
reasonable values of the necessary parameters. 

D. MULTIPLE CORRELATION RESULTS 

To facilitate the overall prediction of the 

heat transfer coefficient as a function of those 

variables considered, the actual shell side coefficient, 

h , was calculated as a product of h. , , and the 

o r ideal 

aforementioned factors: J,, J , J„. Colburn i -factor 
correlations of Table 1 were used for the various flow 
regimes as were correlations for N t and The 

fluid properties and mass flow rate terms were factored 
out of the equation and only those terms involving 
geometric parameters were used for the overall cor- 
relation. These terms were then correlated using the 
aforementioned Fortran subroutines. The correlation 
results are presented in Table 6 for the various flow- 
regimes and tube orientations. It should be noted that 
all input dimensions for the exchanger parameters are 
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in inches and that the correction factors inherent in 
these equations have been accounted for in the cor- 
relation constant. The resulting units for h Q are 
those typically expected of the English system, Btu/HR- 
FT^-°F. The units for fluid properties are as listed 
in the nomenclature. 

Table 6 is a "busy" listing and many graphical 
representations can be drawn plotting the heat transfer 
coefficient h Q against the various parameters. That 
will not be done here but some aspects of Table 6 are 
noteworthy. 

Most exchangers are designed for operation in 
turbulent flow and for this regime it can be seen that 
h Q is nearly inversely proportional to shell diameter, 
the exponent a^ being close to -1.0. This is a striking 
result, in that for desired increases in h , changes in 
shell diameter maybe the most rapid means of attaining 
the desired value. 

The dependence of h Q on baffle cut ratio, 

& c /D^, is independent of flow but dependent on tube 
orientation, and h Q varies , nearly as the square-root 
of l /D.. Pitch- to-diameter ratio, P/d , has little 
effect on square tube orientations, but can cause 
significant changes in triangular and rotated square 
patterns. The dependence of h Q on tube diameter and 
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baffle spacing ratio is straightforward and needs no 
further explanation. 

From Table 6 it is quite easy to predict the 
heat transfer coefficient given the exchanger geometric 
parameters and fluid properties. Only a quick Reynolds 
number calculation is necessary beforehand to determine 
the applicable flow regime. Changes in h Q are pre- 
dicted simply by calculating the ratio of the parameter 
varied raised to the appropriate exponent. This is 
quickly accomplished if only one parameter is varied 
at a given time. For large changes one must be cautious 
that the flow regime remains as desired, i.e., turbulent 
or laminar/transition. 
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TABLE 6 Correlation Results for Overall Shell Side Heat Transfer 
Coefficient, h 
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III. CORRECTIONS AFFECTING PRESSURE DROP 

A. BAFFLE LEAKAGE FACTOR -- R-g 

The pressure drop correction factor for baffle 
leakage, R^, accounts for pressure drop changes in much 
the same fashion that accounts for changes in heat 
transfer. It relates the leakage through the baffle due 
to tube-to-baf f le and baffle- to-shell clearances to the 
geometric parameters of the heat exchanger. Equations 
A-2, 4, 6, 7, 8, 14, 15, 20 and 21 are all used to 
calculate R . Dependence on tube orientation is accounted 
for in equations A-4 and A-6 by the use of the correla- 
tions for N t in Equations (1) and (2) and the tube pitch 
relations of Table 3. 

The non-linear correlation results are presented 
in Table 7 and a representative graph as Figure 11. 

Figures 10 and 11 show a striking resemblance of form, 
which is to be expected. One can logically assume that 
an increase in heat transfer coefficient due to less 
leakage through the baffle will also result in increased 
pressure drops. 



-46- 



R * ■ w 


a l a ? 

(P/d 0 ) C<t c 


a 3 

/D ± ) 3 


a 

<v D i 5 


^ c 

CD ± ) 5 




Layout 


a 0 a l 


a 2 


a 3 


a 4 


a 5 


Square 


.2985 .2217 


1.315 


. 1754 


.4674 


. 1584 


Rotated 

Square 


.328 .1798 


1.147 


. 1475 


. 3954 


. 157 


Triangular 


.3702 .1682 


1.0696 


. 1318 


. 3496 


. 1365 



TABLE 7 Correlation Results for Baffle Leakage Factor 
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FIGURE 11 Baffle Leakage Factor vs P/d Q 
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B. BUNDLE BYPASS FACTOR -- R b 

The bundle bypass factor, R^ , accounts for pressure 
drop corrections in the same manner that accounts for 
heat transfer changes. From Equation A- 22, it can be 
seen to relate that correction to F sb p> the fractional 
crossflow area available for bundle bypass. 

Table 8 presents the correlation results for R^ . 
Similar to is the weak dependence on tube diameter 
and the strong dependence on pitch- to- diameter ratio. 

This comes as no surprise since increases in P/d Q "open 
up” area for flow through the crossflow bank, thereby 
decreasing the fraction of bypass area. The increase 
in R b with increasing is less than that of P/d Q , 
since the shell-to-outer-tube-limit clearance (Equation 3) 
also increases with D^. 

Figure 12 is a graphical representation of R^ 
vs. P/d Q at various D^, in the same manner as Figure 8. 

It should be noted that over the range of d Q considered 
(0.5 < d Q < 1.5), the value of R^ will decrease about 
5% with increasing d Q , a much larger percentage decrease 
than can be predicted for (and subsequently, heat 
transfer) . 
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FIGURE 12 Bundle Bypass Factor, vs P/d Q 
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C. CROSSFLOW AREA PRESSURE DROP 

Equation A-23 gives the ideal pressure drop for 
flow across the tube bank between baffles. When 
multiplied by R^ and R^, the effects of bundle bypass 
and baffle leakage are accounted for. This quantity 
has been correlated in the same manner as the previous 
values. The friction factor correlations of Table 2 
were used along with other corrections for tube 
orientation. Table 9 is the result of these overall 
correlations for the two flow regimes and various 
orientations. Note that the exponent x^ of the basic 
equation depends on flow and orientation as a result 
of the friction factor correlation. Hence, the 
pressure drop can not be expressed only as a function 
of exchanger parameters. From Table 9 it is obvious 
that all the geometric parameters inversely affect 
pressure drop, a result that may have been intuitively 
obvious, but the extent of which can now be reasonably 
quantified. Essentially, the increase in any of the 
five basic parameters allows for more area for fluid 
flow with less restriction; hence, lower pressure 
drops . 
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TABLE 9 Correlation Results for Crossflow Area Pressure Drop 
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D. WINDOW AREA PRESSURE DROP 

Equation A-23 gives the ideal pressure drop for 
the baffle window section of the heat exchanger. When 
multiplied by R^, the effect of baffle leakage is taken 
into account. Correlation results are presented in 
Table 10 for differeing tube orientations. Unlike 
crossflow pressure drop, this quantity is independent 
of flow regime. 

These results are very similar to the crossflow 
pressure drops in that pressure drop is inversely 
proportional to all variables. However, in the window 
area the relative importance of the parameters is slightly 
different. For the window area, shell diameter is the 
predominant influencing factor instead of pitch- to-diameter 
ratio. The relative influence of baffle cut and a baffle 
spacing ratios is nearly the same for both area pressure 
drops. Tube diameter impact is relatively the smallest 
for the turbulent crossflow area pressure drop. 
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IV. PERFORMANCE AND SIZING OF HEAT EXCHANGERS 

A. SENSITIVITY ANALYSIS OF HEAT TRANSFER PER UNIT 

PRESSURE DROP 

For any given application, when determining the 
"optimal" size of a heat exchanger, some feel for the 
overall effect of parameter variations on heat exchanger 
performance is needed in order to judge the direction 
in which the design variables should be altered to arrive 
at a reasonable result. The performance of a heat 
exchanger, or the criterion against which an exchanger- 
is to be considered optimal, can vary widely depending 
on the depth of the analysis that is being performed. 

In most practical applications the ultimate criterion 
is cost. In this analysis, cost will not be considered 
but a factor which can influence cost to some extent 
will be analyzed. 

From the corrleation results for heat transfer 
coefficient, Table 6, and the basic heat balance: 

Q = AU AT. (8) 

o £m 

consider the case where U Q is primarily dependent on h Q . 
Equation (8) then can be written as: 
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Q 



AT 



£m 



W K. 
s 1 



Ah o 



where is defined in Table 6. But A = ird N jV 12 (ft“) , 

and the results from Equation (1) (or (2)) can be used 

for N. to express A in terms of d , D- , P/d . Define: 
t 0 10 



as 



Vi 



which yields: 



V L - 



tt d N h 
0 t 0 

12 



(9) 



This quantity is a non-dimensional expression for heat 
transfer per unit length of the heat exchanger. From 
the correlation results for crossflow and window area 
pressure drops, (Tables 9 and 10), we may form, 
similarly: 



r /L = 



N. 



W /a Q 
S /c, C M S 



(AP 



b,i 



AP 



w 



,i) /L 



( 10 ) 



This is a simplified expression for the total shell-side 
pressure drop per unit length of heat exchanger. Then 
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F h /r p can be expressed as a quotient of those correlations 

already derived. It can be seen from all the correlation 

results that r^/T is not only dependent on the geometry 

of the exchanger but flow rate also, specifically, 

W s /y s . For this analysis, we will consider only turbulent 

flow across square-oriented tube bundles. The applicable 

correlation results were selected, the ratio r,/r 

h P 

formed and that ratio was graphed as a function of the 

exchanger parameters. It should be noted that a Reynolds 
4 

number of 10 was used, resulting in a value of 

W /y = 2.026 x 10 5 at the mid-range of the exchanger 

geometry variables (d = 1.0, P/d Q = 1.5, £ S /D^ = -6, 

D\ = 35). Figures 13, 14, 15, 16 and 17 are various graphs 

of T, /T . From Figures 13-17 it is obvious that shell 
h p 

diameter is the most influential parameter affecting 
heat exchanger performance. Figure 13 graphically 
illustrates this influence for several values of baffle 
spacing ratio, £ s /D^. This is not to say that we should 
always optimize an exchanger at large values of D^. 

Although r h /T p increases, this graph does not allow 
for any constraints to be considered, e.g. maximum shell- 
side pressure drop. However, we may conclude that where 
large amounts of heat are to be transferred with low 
shell-side pressure drops as a requirement, "large” shell 
diameters will be required. 
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Figure 14 illustrates the influence of baffle 
spacing on F^/r^ at several values of D^. This parameter 
also favorably affects the heat exchanger performance 
but certainly not to the extent of shell diameter. Again, 



the realization of the absence of a pressure drop con- 
straint must be accounted for. 

Figure 15 illustrates the relationship between 
and tube diameter. This is the only parameter (of 
those considered) which adversely affects r^/F^. The 
most basic explanation for this effect is that increases 
in tube diameter will more quickly decrease heat transfer 
than they will pressure drop. 

Figure 16 and 17 relate the effects of baffle 

cut ratio, l /D . and pitch- to-diameter ratio, P/d , to 

T, /T . It is obvious that at the conditions analyzed 
h p 

(turbulejit flow, square tube orientation) , neither para- 
meter has any significant effect on heat exchanger 
performance . 

The five figures shown can be of assistance 
when determining exchanger size for a given "job". The 
qualitative effect (or lack thereof) of the five para- 
meters will help guard against "false starts" in any 
iterative sizing scheme. 
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FIGURE 15 I\/r vs D. 

n p l 
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FIGURE 14 r h / r p vs V D i 
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FIGURE is r,/r vs d 
n p o 
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B. PROPOSED SIZING ROUTINE' FOR CASE OF SHELL -SIDE 

PROPERTIES ONLY 

For the case where the shell-side heat transfer 
coefficient is the controlling factor (the tube side 
coefficient has little impact on overall U ) , it is 
desirable to determine the exchanger geometry for a 
given "job": total Q transferred ; AP specified; W , 

AT g (shell-side) , and fluid properties all given. We 
may desire to "optimize" the heat exchanger on minimum 
volume where volume is simply expressed as: 



4 



From the basic heat balance equation, 

tt d N L 

Q = AU AT, , U = h , A = — - — (11) 

x o £m oo ^2 



yields 



tt d N. L 

Q = AT. — - — h 

2,m ^ j o 



( 12 ) 



With knowledge of the fluid properties we may form the 
quantity 



IV = 



AT. K-, W 
Jim Is 



(13) 
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where is specified in Table 6, correlation results 
for shell-side coefficient. Then 



V L 



* d Q N t L 
12 



h 



o 



(14) 



We must now begin to specify the geometry of the 
exchanger. First, pick d Q , P/d Q , as the beginning 
points. From equation A- 12 for Reynolds number:- 



R 



12d W 
_ o s 

s u s m 
s m 



(15) 



W 3 . 3. a. 

- — =0 (V 1 ( p/d 0> 2 'W 3 tv" 4 

^s 



where 



c = 62.224, a, = 1.037, a- = - 1.915, 
o 1 2 



a- = 1 . 0 , a. = - 1.935 
o ’ 4 

we may either determine J& s /D^ (if we desire to operate 
at a certain R„ ) or pick l /D. and determine the flow 
regime (laminar/transition or turbulent) . 

Once this condition has been finalized, the 
proper correlation for h Q is selected from Table 6, 
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knowing the selected tube orientation; The last para- 
meter to be determined is £ C /D^. As Bell has stated, 
for "well-designed" exchangers, J c should be about 1.0, 
which, from Figure 6, yields a value of £ C /D^ of about 
0.25. The initial value of h Q may now be calculated 
from the correlation results of Table 6 and exchanger 
length L determined from Equation 10, using the proper 
correlation for tube orientation. Volume may then 
be easily calculated. 

Pressure drop calculations can be made from the 

correlation results of Tables 10 and 11, knowing the 

exchanger length L. This must then be checked against 

the allowed AP . If AP . < AP , we may 

s calc allowed 

decrease the size of the exchanger to utilize more of 
the allowed AP g . This is most easily accomplished by 
decreasing and recalculating h Q , L, AP g and V, to 
arrive at some minimum volume. 

If AP,. , > A P c ,, ,, the procedure is 

s calc s allowed ^ 

reversed, i.e., increase and recalculate the 
necessary quantities until a convergence within some 
desired bound is obtained. As in the case in any 
sizing procedure, some notes of caution and guidelines 
are in order. When increasing shell diameter (or any 
of the other parameters) , caution must be taken that the 



-67- 



flow regime remains as desired. It is also possible to 
"fine tune" the exchanger by altering the other geometric 
values. If the difference APs ca i c * AP s allowec i small, 
we may desire to change the baffle cut or baffle spacing 
ratios. Small changes in these parameters will also 
result in changes to the shell side pressure drop and 
heat transfer coefficient, the effect of which can be 
examined by the exponents of the applicable correlation 
result. We may also use Bell's insight into the heat 
transfer corrections to assure that the heat exchanger 
is reasonably well designed. For the heat transfer • 
correction factors , J c> Bell gives the following 

ranges as good guidelines: 

J * 1.0 

c 

0.7 < J, < 0.9 

D 

0.7 < < 0.8 

From Figure 9, 10, and 11, we can pick reasonable values 
of and £ s /D^ as starting points for the exchanger 

geometry and, should changes be made, ensure that we 
remain in a normally expected range for those values. 

The use of those figures as aids to sizing should prevent 
haphazard variations and point to a reasonable solution 
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FIGURE 18 Simplified Logic Sequence for 
Proposed Sizing Routine 
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in an expedient manner. Figure 13 gives a simplified 
flow chart of the logical sequence for this proposed 
sizing routine. 

C. EXAMPLE OF PROPOSED SIZING ROUTINE WITH 

' COMPART SON TO THE DELAWARE METHOD 

The preceding sections of this analysis have dealt 
with the actual results of the non-linear regression used 
to simplify the heat transfer coefficient and pressure 
drop calculations. In this section, a hypothetical 
example for determining exchanger size is given to show 
the actual mechanics of the proposed sizing routine and 
to compare those results with the calculations of the 
Delaware Method. This comparison is done for the case 
where only the shell-side heat transfer coefficient and 
pressure drop are of concern to the designer. 

The "job" considered in this example is as 

follows : 

1) shell-side mass flow rate of water: 

W s = 250,000 lbm/hr 

2) heat water from 50°F to 90°F 

3) maximum shell-side pressure drop: 

AP S <_ 8.0 psi 
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The following fluid properties are assumed: 



1) c « 1.0 btu/lbm - °F 

2) y , = 1.65 lbm/hr-ft (at an assumed 

S y W 

T wall = 100 °^ 

3) y s = 2.445 lbm/hr-ft (at T ave = 70°F) 

4) P rs = 7.05 

5) p s = 62.075 lbm/ft 3 

We desire to operate in the turbulent flow regime 
and must check the Reynolds number for our exchanger. 
Picking tube diameter d Q = 0.75 in, pitch- to-diameter 
ratio P/d Q = 1.33, baffle spacing ratio £ s /D^ = 0.75 
and shell diameter = 25.0 in, a check of Reynolds 
number using Equation (15) yields Re s = 7229, which 
ensures turbulence. Using a square tube orientation, 
baffle cut ratio i, /D^ = 0.25, we select the applicable 
correlation result from Table 6: 



h 



fw ' “ ‘ ^ 



- .212 



.015 



o 



(4.07) (d Q ) 




-.4492 



-.376 



-1.032 






whe re 



0.14 



-2/3 y 



s 



K, = c (P ) 

1 S 
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Using the values given results in a heat transfer 
coefficient of: 



h Q = 463.5 Btu/hr-ft 3 -°F 



From Q = AU AT = Ah AT' , where 
o Jin o x.n 



2.17526 



f d N t L D. 

A = — - — and N = 0.3559 — 7 

12 r d Z 



The length L of the heat exchanger is calculated: 



L = 6.99 f t , 
then volume 

V = 23.83 ft 3 . 



Using the correlation results for crossflow and window 
area pressure drops from Tables 9 and 10 respectively, 
yields the total shell-side pressure drop: 

AP = .29 psi 
s 

As a comparison, utilizing the Delaware Method (Reference 
1) , the following is presented: 
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Proposed Method 


Delaware Method 


h Q (Btu/hr- ft' -°F ) 


46 3.5 


454.8 


L (ft) 


6.99 


7.05 


V (ft 3 ) 


25.83 


24.03 


AP (psi) 


0.29 


0.25 



Although the results are quite close, we are far from 
"optimized" in using the allowed shell-side pressure drop. 
We may decrease shell diameter to most easily 
accomplished this. For the next interation, pick 
D. = 15.25 in. This "odd" size is picked only because 
Bell's tables in Reference (1) list all the applicable 
parameters for this size (D Qt ^, , etc.) and inter- 
polation between those listed values would not result in 
a meaningful comparison. 

A recalculation of the heat exchanger using this 
new value of yields the following result: 





Proposed Method 


Delaware Method 


h o 


785.16 


789.0 


L 


12.09 


11.69 


V 


15.34 


14.83 


AP 


3.62 


3.61 
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These results are as agreeable as those of the 
first iteration, yet the exchanger may be further de- 



creased in size. 


A third iteration. 


using = 13.25 in 


yields : 








Proposed Method 


Delaware Method 


h 

0 


907.8 


917.5 


L 


14.20 


14.61 


V 


15 . 6 


13.99 


AP C 


7.53 


7.27 



These results show that we have come close to a 
preliminary size for the heat exchanger that meets the 
required "job". Some comments concerning the previous 
calculations are noteworthy. 

In decreasing shell diameter in these iterations, 
the recalculation of h Q by the proposed method is 
accomplished simply by multiplying the previous value 
by the ratio of shell diameters raised to the applicable 
correlation exponent. Pressure drop calculations are 
estimated in much the same manner, except that the 
crossflow and window area terms must be done separately. 
The only "lengthy" calculation involves that of exchanger 
length, L. 
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For the sizing determination by the Delaware 
Method, no less than twenty-one (21) separate calculations 
and seven (7) graphical readings must be done for each 
iteration. For this proposed method only five (5) 
calculations were necessary: .heat transfer coefficient 

(Table 6) , overall length (Equation 12) , pressure drops 
for crossflow and window sections (Tables 9 and 10) , 
and total pressure drop (Equation A-25). 

This quick, handy method is, in the opinion of 
this author, much more preferable to the original method 
for preliminary heat exchanger sizing and "back-of- the- 
envelope" calculations. It is realized that no simplified 
heat exahanger model can replace the more sophisticated 
computer optimization methods that have been and are 
being developed. However, in a relatively short time, 
the designer can achieve a realistic estimate of the 
heat exchanger performance and size, and can use these 
results as possible inputs to the more elaborate routines. 
The closer an initial estimate is to the actual optimum, 
the more time and money can be saved in the more 
sophisticated (and expensive) methods. Minimization 
of cost will almost always be a major goal. 
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APPENDIX A 



Basic equations used for heat transfer coefficient 
and pressure drop calculations. 



(1) N = 



2jl c 

D i ci - Tjf: 



( 2 ) 



F = — 

C TT 



IT + 2 



2 cos 



-1 




sm 



cos 



, D. - 21 
-1 i c 



D 



otA 



[ ot l J 



(3) N 



0.8 l 



cw 



(4) 



S = l 
m s 



^i ^otJi + 



fD . - d } 
ot& o 



n 



(P - d 0 5 



(5) 


F sbp 


■ CCDj • 


• D otJ.^ 


V s . 


(6) 


S tb = 


* d o 4 tb N t < 1/2 > d 








21 1 




(7) 


3 = 2 


COS 


1 - C 










D i> 
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C8) S sb = */2 Dj 6 sb Cl - 3/2] 





D. 2 




21 } 


— 


(9) 


S = — 

W S 4 


3 

2 


1 - — £ 
D. 

l J 


sin 3/2 



C 10 > S wt “ f (1 ' F c>* d o 2 
C n > S w ' S wg * S wt 



d W 

(12) R = - 2 —± . 12 

e s y S 







fw 1 


-2/3 f 


(13) 


^ideal ^"i c s 


s 

s 

m J 


CP rs ) 


(14) 


S , + S«, 

sb tb 






T lm n 

m 







s 



> 






r 



s 






. 14 



(15) 
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(16) 



(17) 



(13) 



(19) 



( 20 ) 



( 21 ) 



( 22 ) 



J = 0.55 + 0.72 F 
c c 



J A - 0.44 (1 - r s ) + Cl 



2 . 2r 



0.44(1 - r s )]e 



Jim 








o ideal c l b 



p = - 0.15 (1 + r ) + 0.8 



R * = 



[-1.53(1 + r s )rP m ] 



R, 








AP 



b , i 



4 ^ W c 

is c 

2 p g S 2 
H s & c m 



[ y 



0.14 



s ,w 



(23) 



2 
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(24) 



( 25 ) 



AP • = 
w,i 



W/ (2 + 0.6 N ) 
s cw 

2 g p S S 
6 c K s m w 



AP = C(N, + 1) AP, - R, + N, AP . ]R 0 
s '■b ' b,i b b w , i £ 
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( 2 ) 

(3) 

(4) 

( 5 ) 

( 6 ) 

(7) 

( 8 ) 

(9) 

( 10 ) 
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